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We investigate within a self-consistent theory the molecular instabilities arising in the normal 
state of a homogeneous degenerate Fermi gas, covering the whole BEC-BCS crossover. These are the 
standard instability for molecular formation, the BCS instability which corresponds to the formation 
of Cooper pairs and the related Bose-Einstein instability. These instabilities manifest themselves in 
the properties of the particle-particle vertex, which we calculate in a ladder approximation. To find 
the critical temperatures corresponding to these various instabilities, we handle the properties of 
the interacting Fermi gas on the same footing as the instabilities by making use of the same vertex. 
This approximate treatment is shown to be quite satisfactory in a number of limiting situations 
where it agrees with known exact results. The results for the BCS critical temperature and for the 
BE condensation are found to be in fair agreement with earlier results. The threshold for formation 
of molecules at rest undergoes a sizeable shift toward the BEC side, due to quantum effects arising 
from the presence of the degenerate Fermi gas. This should make its experimental observation fairly 
easy. This shift remains important at least up to temperatures comparable to the Fermi energy of 
the gas. 

I. INTRODUCTION 

Experimental progress in the field of the BEC-BCS crossover in ultracold fermionic gases has been going on recently 
at a very fast pace [1] . In particular after the observation of the Bose-Einstein condensation of diatomic molecules [2] 
both in 6 Li and in 40 K, evidence for superfluidity on the BCS side of the crossover has been provided by the study of 
collective modes by Bartenstein et al [3], which shows a strong attenuation peak, whose likely interpretation is pair 
breaking in the vicinity of the BCS transition. Further evidence have come from "projection" technique experiments 
[4] where a fast sweep from the BCS to BEC side allows to infer the existence of superfluidity on the BCS side. Very 
recently the observation of vortices in these Fermi gases on the BCS side [5] has given a much clearer evidence for 
superfluidity. 

One major interest of the study of the BEC-BCS crossover is to obtain a clear and precise picture of the way in 
which Cooper pairs go progressively into diatomic molecules when the strength of the attractive interaction between 
fermionic atoms is progressively increased, and in particular to explore if the BCS formalism gives a proper description 
of this evolution. The occurence of this crossover within the BCS formalism and its interest has been put forward 
by the works of Leggett [6], and Nozieres and Schmitt-Rink [7], and further investigated by Sa de Melo et al [8], 
although the fact that the BCS formalism gives also a correct description of the condensation of molecules in the 
strong coupling limit of dilute molecules was known much earlier [9,10]. 

Apart from its fundamental physical interest this question is also highly relevant for high T c superconductors since 
in this case Cooper pairs are known to be quite small and much closer to molecules than in standard superconductors. 
The situation in high T c superconductors would accordingly be close to a Bose condensation of molecules, which would 
provide a natural explanation for the high value of the critical temperature. The natural simple theoretical framework 
for this solid state problem is the attractive Hubbard model and the self-consistent T-matrix approximation [11,12] 
used to obtain a qualitative analytical understanding of the Hubbard model corresponds exactly to the framework we 
will use in the present paper to handle the case of fermionic atomic gases. 

In practice the BEC-BCS crossover is realized experimentally by going through a Feshbach resonance, by varying 
a static and homogeneous magnetic field applied to the atomic gas. Due to the very low temperature only s-wave 
scattering occurs in the gas and one needs atoms belonging to two different hyperfine states in order to obtain a 
non-zero scattering, which is forbidden by Pauli principle between atoms belonging to the same hyperfine state. The 
scattering length a varies accross the Feshbach resonance, starting with fairly small negative values at high field, which 
corresponds to a weakly attractive effective interaction between atoms belonging to two different hyperfine states. 
When the field is lowered a goes to highly negative values and diverges right at the resonance, where it jumps an 
infinite positive value. When the field is further lowered the positive scattering length a decreases down to fairly small 
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values. On the low field side where a > 0, molecules exist made of two fermionic atoms belonging to two different 
hyperfine states, while such a bound state does not exist on the a < side. 

An experimental surprise has come from the fact that this singular behaviour, which occurs for two atoms in vacuum 
or for dilute Fermi gases, disappears when one goes to dense gases. For example the energy of the gas, measured 
in expansion experiments [13], does not display any singularity and is perfectly smooth when one goes through the 
resonance, while the scattering length a displays a singularity at the same location. This has been explained [14] by 
the effect of the dense Fermi gas, which introduces another length scale, namely the Fermi wavelength, related to 
the Fermi wavevector kp, defined from the density n of atoms belonging to a given hyperfine state (we assume that 
the density of the two hyperfine states are equal) n = k F /(6Tr 2 ). The associated energy scale is the Fermi energy 
Ep = k F /{2m). The presence of the dense gas washes out the Feshbach resonance. 

This effect is quite reasonable when one recalls that the very existence of the Cooper pairs, in the case of a 
very weak interaction, is due to the existence of the Fermi sea without which a bound state could not form. In 
other words the existence of the dense fermionic gas, or equivalently the Fermi sea, affects the formation of bound 
states, that is of molecules. In particular it has been shown explicitely that [15], in the normal state, the threshold 
for the formation of molecules is affected by the presence of the dense gas. Instead of having, as in vacuum, all 
the molecular bound states appearing (with zero binding energy) at the same magnetic field corresponding to the 
Feshbach resonance, characterized by a -1 = 0, the threshold for the appearance of these bound states depends now 
on the total momentum of the molecule. Qualitatively this effect is easily understood. Indeed, seen from a k-space 
point of view, the formation of the molecular bound state requires the partial occupation of plane-wave states (the 
amplitude for the occupation probability is the Fourier transform of the molecular wave function). However in the 
presence of a dense gas, some of these states are already occupied and Pauli exclusion makes them unavailable for 
building up the molecular wave-function. As a result, with less plane-wave states available, the formation of the 
molecular bound state becomes more difficult and the threshold is naturally pushed toward stronger interactions, that 
is toward positive scattering lengths a. 

It is also clear that this effect depends on the total momentum K of the molecule. If this momentum is very high, the 
molecular formation is unaffected since the plane-waves required to make up the wave function are near K, and they 
are all essentially empty. In the case the threshold is at the same location as in vacuum, that is a -1 = 0, independent 
of temperature. On the other hand for a molecule with K = the required plane-waves have small wavevectors and 
these states are partially occupied. This is the case where the presence of the dense gas is the most strongly felt and 
the shift of the threshold is strongest, and this is the case we will more specifically consider. Naturally in all cases the 
effect depends on temperature since, when temperature is raised, the occupation of plane-wave states gets smaller and 
accordingly the effect on molecular formation is reduced. Ultimately at quite high temperature, the effect disappears 
because all the plane-wave states have very low occupation probability and we go back to the classical gas situation 
where the threshold is at a -1 = for any molecular state. 

Naturally the experimental observation of this effect would be quite interesting, not the least because it goes against 
the simple intuition one gets from classical gas physics. This implies to have reasonable evaluations of the domain 
of physical parameters where it occurs. Unfortunately the evaluation made in Ref. [15] was very rough since the 
parameters of the Fermi sea were taken as those of a non interacting Fermi gas. This is naturally quite inconsistent 
since the appearance of molecules is due to interaction, which are in particular expected to have strong effects in the 
vicinity of the Feshbach resonance. This is even more so when one goes toward the BEC side, since one ends up with 
a system which is physically described as a dilute gas of molecular bosons, which has naturally little to do with a free 
Fermi sea. The initial purpose of the present paper is to proceed to a coherent calculation and describe the Fermi 
sea, taking into account self-consistently the interaction responsible for molecular formation. 

More specifically the location of the molecular threshold has been obtained in Ref. [15] by writing that the full 
vertex for atomic scattering has a pole at zero energy ui = 0. This vertex has been obtained by a ladder approximation, 
which is a natural extension of the exact treatment for two atoms in vacuum. Here the existence of the dense gas has 
been taken into account by making use of the free fermion propagator for a given chemical potential /i, instead of the 
fermionic propagator in vacuum. Now if we want to have consistently the properties of the fermionic gas, we should 
make use of this same vertex within the ladder approximation as a starting point. This is what is done in the present 
paper. From the vertex, we will obtain the fermionic self-energy, and then the fermionic density n as a function of 
chemical potential \i and temperature T . In this way we will have in a consistent manner what is essentially the 
equation of state of the fermionic gas in the presence of interactions. 

The appearance of the molecular bound state is not the only feature which is signaled by a pole in the vertex, the 
BCS transition appears also as a pole for an energy per particle equal to the chemical potential, so uo = 2/i. This 
corresponds physically to the appearance of Cooper pairs, which arc quite analogous to molecules. Hence it is a natural 
development to consider also the location of the BCS transition within our approximation. It is obviously of high 
interest to go toward the BEC side of the crossover and we will accordingly also determine the critical temperature 
for the Bose condensation of molecules, once we are beyond the threshold for their formation. All these features of the 
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phase diagram, namely appearance of molecular bound states, BCS transition and Bose condensation of molecules, 
have the common property to be instabilities of the normal state and are revealed by the study of the same vertex. 

As it happens the above approximation, we have been lead naturally to consider, has been already quite often 
considered in the literature in other contexts, dealing with interacting Fermi systems. It is often named the self- 
consistent T-matrix approximation and, as we mentionned above, it has been in particular used to handle the Hubbard 
model (see [11,12] and references therein). In the context of ultracold fermionic gases, it has been already studied 
and used, in particular in the supcrfluid phases (see Ref. [16-18] and references therein). More specifically Pieri and 
Strinati [16] obtain their equations as the result of a specific regularization for the well-known ultraviolet divergence 
which appears in the theory for a contact interaction. In their way the omission of diagrams other than ladder 
diagrams appears as an exact result. In our case we do not restrict ourselves to a contact interaction and keeping 
only ladder diagrams is just an approximation (physically quite reasonable), considered as to be possibly improved. 
Nevertheless we end up with the same equations, as can be more clearly seen by taking the normal state limit of the 
equations of Perali, Pieri, Pisani and Strinati [17]. 

On the other hand we will be concerned only by the normal state properties in the present paper, with more 
specifically in mind the study of the threshold for the molecular bound state, which might be relevant for quite recent 
experimental work [5] as we will see. On the other hand their work [17,18] concentrates on the superfluid phase. Also 
they have focused on trapped gases, while we will deal only with the homogeneous gas. Actually the only overlap is 
at the level of the critical temperature of the superfluid phase, which they have also calculated for the homogeneous 
system [17]. However our specific handling of the equations are different from theirs, each one having its own interest. 
More details and comparisons will be given in the course of our paper. Another specificity of our paper is that, in order 
to check the quality of our approximation, we will compare it in various non-trivial limiting cases with exact results 
already known in the literature. In all cases we will find a perfect agreement which establishes the self-consistent 
T-matrix approach as a physically very reasonable approximation. To our knowledge this specific comparison has not 
yet been done. Yet it is very important to assert precisely the quality of our description of the normal state. This is 
obviously quite necessary if we want to have a proper description of the superfluid state. This point has already been 
stressed quite a number of times in the context of high T c superconductivity, where it is often asserted that the normal 
state properties are even less properly understood than the superconducting ones. We stress that our handling of the 
equations will be made without any approximations, with either analytical or numerical treatment. In summary the 
present paper is complementary to the work of Pieri, Strinati and coworkers. In particular it shares the same spirit of 
trying to describe the physics of these ultracold gases starting from a single coherent fermionic picture for the whole 
crossover, in contrast to more phenomenologically oriented approaches where a bosonic field describing the molecular 
states is introduced immediately at the level of the Hamiltonian. 

The organization of our paper is as follows. In the next section we introduce our basic equations to calculate 
the particle number in terms of the chemical potential and the temperature and we give explicitely the expressions 
for the self-energy and the vertex we will use throughout the paper. Then in the following section we consider the 
high temperature range and show that in this regime our approximation reduces to the exact virial expansion of 
Beth and Uhlenbeck [19,20] for the quantum gas. The interest of this limiting case has been emphasized recently by 
Pitaevskii and Stringari [21] who showed that it does not have any singular behaviour at the unitarity limit when 
the Feshbach resonance a^ 1 = is crossed, despite the appearance of molecular bound states. Then we show that, 
in the large momentum regime, at zero temperature, our approximation reduces to the known exact asymptotic 
result of Belyakov [22] . It is of particular experimental interest to emphasize this limiting case since very often the 
high momentum tail of the particle distribution is analyzed experimentally in order to extract the temperature, by 
comparison with the ideal case distribution, thereby ignoring the effects of interactions on this tail. We analyze also 
the dilute regime on the BCS side and we show how our approximation can be corrected to obtain perfect agreement 
with Galitskii's result [23]. Finally we show explicitely how the molecular bound state arises from the vertex and give 
the corresponding contribution to the self-energy. We analyze then the dilute limit on the BEC side of the crossover 
and show that one ends up, in the normal state, with the expected physical situation with the proper Bose distribution 
for molecules and the expected wave function for this nearly resonating situation. Finally we discuss in this regime 
another temperature, which is quite important in the context of high T c superconductivity, namely the pseudogap 
temperature T», which corresponds to the smooth crossover between the low temperature domain, where molecules 
(i.e. preformed pairs) dominate and the high temperature range where one has only free fermions. Then, gathering 
the different contributions to the self-energy, we calculate the momentum distribution for the atoms and total particle 
number. We emphasize in particular that a separation between particles corresponding to free atoms and particles 
belonging to a molecule has a quite restricted range of validity, although this analysis is quite often found in the 
literature in phenomenological analysis. Finally we end up by displaying and discussing our phase diagram for the 
molecular, the BCS and the Bose-Einstein instabilities. 
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II. SELF ENERGY AND PARTICLE NUMBER 



since 
In our 



The particle number n in a single hyperfine state is obtained from the temperature Green's function G(k, ioj n ), 
where iv n = (2n + l)nT (with n an integer) is the Matsubara frequency, by (we take Ti = 1 and ks = 1): 

n = T J2 j^ G ^^^ nT w 

where r — > + . It is convenient, for actual calculations, to separate out in this equation the free particle contribution 
Go(k, iu n ) = [iuo n — (ek — where £k = k 2 /2m is the free particle kinetic energy, for which the result no for the 

particle number is known to be: 

„ f°° dk fc2 
n = 4ir / . ., j (2) 

Jo (2tt) 3 exp(%) + 1 

where we have set £k = — fj,. The remaining contribution is given by: 

^ T \p [ dk Y,(k,iuj n ) 

" H ° (27T) 3 [iw n -&-£(fc,iWn)].[iWn-&)] W 

the Green's function is related to the self-energy T,(k,ioj n ) by G(k, iw n ) = [iw„ — £k — S(fc, iw„)] _1 . 
our ladder approximation the self-energy can be written as: 

£(fc,io; n )=:r£ J ^T(K,iw v )G (K--k,iw v -iw n ) (4) 

In contrast with the more general situation [24] , we have here a single term because only atoms belonging to different 
hyperfine states do interact, which forbids an exchange term. Here Y(K.,iuj v ) is the standard vertex [24] in the ladder 
approximation. Note that in this vertex, u>„ = 2-kvT (with v being an integer) is a bosonic Matsubara frequency. 
Since the wavevector dependence of T(K,iu v ) is actually only on K, the angular integration of Go on K can be 
performed in Eq. (4) and gives: 

f , n n , v , ■ \ 27rm iuj m + fi- K 2 _/2m 

/ dn K G K - k,iuj m ) = — — In- ■ (5 

J kK iu) m + pt - K^_/2m 

with K ± = K±k. 

In the following we find it more convenient to work with reduced variables. We take as unit of energy the absolute 
value of the chemical potential |^|. This is the natural choice at low temperature in the degenerate regime, but not 
at high temperature in the classical regime where T is a more convenient scale. Anyhow it is easy to switch in our 
formulae from one scale to another. Another inconvenience with our choice is that we have to keep track of the sign 
s of /Lt, defined by s = = ±1, since it switches sign when one goes from the degenerate to the classical regime. 

We similarly define a wavevector scale ko by fcg/2m = \fj,\. Hence we set u — \/j,\oj and so on for all the frequencies, 
together with k = kok and K — k K. We also introduce the reduced temperature t = T/\fj,\. 

The explicit expression for T(K,iuj u ) is obtained, for example, as in Ref. [14], but keeping a nonzero wavevector: 

r 1 (K,^) = -^- + J ^[Tj2Go(^',iuJ m )G (K--k',iw v -iuj m )-^-} (6) 



where we have introduced the standard coupling constant A = —2k a/n. The angular integration and the Matsubara 
summation can be performed and give, with the change of variable k' = k x: 



T-\K,u) = X- 1 + / dx 



f 

Jo 



xt/K cosh [((a; + A^/2) 2 - s)/2t] 



x 2 + (K/2) 2 -s-u/2 cosh [((a; - K/2) 2 - s)/2t] 



(7) 



where we have introduced a "reduced" vertex defined by f = — mk T / (2tt 2 ) and used reduced variables. 

In principle one could keep working with Matsubara frequencies and perform numerically the discrete frequency 
summation coming in the expression of the self-energy Eq. (4) . This is the procedure chosen for example by Perali et al 
[17] . However the corresponding series turn out to be rather slowly converging, which is numerically unpleasant. Hence 
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we will rather transform here this summation into integrals over frequency in a standard way [25] by introducing the 
Fermi distribution f(uj/T) = l/(exp(w/T) + 1) which has its poles at iu> n , writing the summation as an integral over 
contours encircling these poles and deforming the contours. The contribution coming from the part of the contour at 
infinity is checked to be zero, and in general we are left with two contributions. We leave out for the moment a third 
possible contribution from a pole of T on the real negative frequency axis, corresponding physically to a molecular 
state. This will be taken up in section VI. Both of our contributions turn out to be rapidly convergent for large 
frequencies because of the Fermi distribution. One of them, which we call Er, arises from the cut of T(K, iuo n + w) 
which extends from the branch point <I>b = 2((K /2) 2 — s) — ioj n to oo. The other one E^ comes from the logarithm 
in Eq.(5) and encircles the cut going from w_ = K 2 _ — s to uj + = K? — s, where K± — K ±k. For both contours we 
can rewrite their contributions as an integral along the corresponding cut by introducing the jump across the cut of 
the function to be integrated. We obtain in this way T,(k,iuj n ) = E(fc, iuj n )/\^i\ — E r (/c,wl> n ) + Si(fc,zw„) where: 

S L (A,iw n ) = --= / dKK / dyf(^)f(K,y + iu> n ) (8) 
2k Jo J Q _ t 



£r(Mw„) = --y dKK dy &(f ) In 

27Tfc J Q JK 2 /2-2s 1 



y-Kl 



y - K+ + S - lLU n 



lmT(K,y) 



(9) 



where ImT(K, y) is for Imr(_ftT, y + ie) with e — > + . We have introduced in the expression for Er the Bose distribution 
b(y) = l/(cxp(y) - 1). 



III. HIGH TEMPERATURE RANGE 



It is of interest to consider the high temperature limit of our approximation. We will sec that it reduces to the virial 
expansion of Beth and Uhlenbcck [19,20] which means that it becomes exact in this limit. This is naturally quite a 
satisfactory feature of our approximation. For simplicity we restrict ourselves to the a < range where there is no 
molecular bound state. We have also proceeded to this comparison in the a > case, where the presence of bound 
states modifies the Beth and Uhlenbeck result, and we have found that our approximation reduces also to their result 
in this case. 

In this limit we go to the classical regime where [i — ► — oo with 1/t = \[i\/T — > oo and in the above formulae we 
have s = — 1. In this case the cosh's in Eq.(7) may be replaced by exponentials and the integral on the right-hand side 
reduces to (n/2)^(K/2) 2 + 1 - Q/2. This is just what one would get from the scattering amplitude of two particles in 
vacuum. This is naturally expected since, in this classical regime, one goes in a dilute limit. One can then see that the 
cut contribution E r to the self-energy goes to zero as e~ 2 \v\l T because of the presence of the Bose distribution b(y/t) 
and of the lower bound on the y integration. It is therefore negligible compared to the pole contribution E^ which 
goes as e~\^/ T as we will now see. Indeed in Eq.(8) we may replace the Fermi distribution f(y/t) = l/(exp(j//i) + 1) 
by its classical limit exp(—y/t). Moreover, since t — > 0, y is restricted in Eq. (8) to a vanishingly small range above 
the lower bound w_ = (K — k) 2 + 1, so in T we can replace y by this lower bound. The y integration is then easily 
performed. Similarly we want to have this lower bound as small as possible, in order to pick the dominant contribution 
from the Fermi distribution. This restricts K to a vanishingly small range around k. Hence we can replace K by k, 
except in the exponential, and the remaining integral is again easily performed. Finally it is more convenient to go 
back from our reduced variables to the physical ones, since in this limit the actual energy scale is (2mT) 1 / 2 . This 
leads to: 

-x 2 ^Q/ 9 ImI 1 47mn/m . „. 

£(fc,^„)=- v -T 3 / 2 e-- 1/2 T7 ^^ = = — / !\ i • ^ ( 10 ) 

Actually we have omitted a term (l/2)efc in the square root, which turns out to be negligible since we will naturally 
find that only efe <~ T is relevant in our calculation, while — ► oo. Hence the self-energy depends only on the 
frequency. 

Let us first consider the case of small coupling |A| <C 1, where the square root in the denominator is neg- 
ligible compared to m _1 / 2 a _1 . In this case the self-energy is just a frequency- independent constant E = 
a(2m/7r) 1//2 T 3 / 2 exp(— \fi\/T). Actually this result coincides, as it should, with the mean field expression E = gn a 
with the coupling constant g — Aira/m. Hence this is just as if we had free particles, and had shifted the chemical 
potential from ji to jj, — E. The resulting change Snk in the single hyperfine state particle distribution is, in this 
classical regime: 
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5n k = --e — (11) 

After integration over k, this leads to a change of particle density: 

Sn = - ^ A T 3 e-W (12) 

where we have introduced the thermal de Broglie wavelength At = (27r/(mT)) 1 / 2 . After substitution of the expression 
of S, this is Sn = — 2a A^ 4 cxp(— 2|/i|/T), which is just the result of the virial expansion for our case. 

Let us now come back to the general case where we have to take into account the frequency dependence of the 
self-energy. For simplicity we keep assuming that a < 0, so we have no bound state. Since the self-energy is small in 
this regime, the integrand in Eq. (3) reduces to iu> n )/{iuj n — £k) 2 - The integration over k is done first and gives 
2i7r(m/2) 3 / 2 / (iu> n — l/^l) 1 / 2 where the determination of the complex square root is with a positive imaginary part. On 
the other hand the square root (|/i| — iujn) 1 / 2 coming in Eq. (10) for the self-energy has the sign of its imaginary part 
opposite to the sign of w„, as it can be seen from the starting expression Eq. (7). Hence this square root has a cut 
on the real positive axis, starting from uj — 

The remaining summation over Matsubara frequencies is transformed in the standard way [25] into an integral over 
frequency on a contour encircling the poles of the Fermi distribution f(u/T) = l/(cxp(w/T) + 1). This contour is 
then deformed into a contour C which goes in a clockwise way around the cut [|^|,oo[ of the square root [26] , which 
gives: 

Sn = -2^{j) J c ^v^TT (uj - |mI) 1/2 (13) 

Taking into account that the determination of (uj — l/il) 1 ^ 2 in the denominator changes sign when one crosses the cut 
[|^|, oo [, this integral is twice the integral of its real part on the cut. Making the change of variable uj = + p 2 /m 
and using Eq. (10) for £(w), we obtain finally in this classical limit \n\/T — > oo: 

5n = - l -a (l^yV'H/T f°° dp f^- (14) 

which is just the Bcth-Uhlenbeck result for our case. In the unitarity limit where mTa 2 — > oo the integral reduces 
to 7r/(2|a|) which gives Sn = (l/2)(mT/7r) 3 / 2 e~ 2 \^/ T . Naturally we have just written here the contribution coming 
from the interaction. The statistical correction [19,20] in the virial expansion will come out if we expand the Fermi 
distribution in Eq. (2). 



IV. LARGE MOMENTUM BEHAVIOUR 



Another limiting situation which is of interest to consider is the large momentum regime. A perturbative calculation 
in this limit has been done a long time ago by Belyakov [22]. Since his calculation was performed at T — 0, we will 
also restrict ourselves to this case. Just as in the above section it is convenient to transform the summation over 
Matsubara frequencies in: 

n k =T^G(k,iw n ) e l ^ T (15) 

n 

into an integral over frequency on a contour encircling the poles of the Fermi distribution and to deform this contour 
into the contour C — C\ + C2 with C\ =]ioo — e, — ioo — e[ and C2 =] — ico + e, zoo + e[ with e — ► + . At T — the 
second part C2 of this contour does not contribute, and we close the first part C\ on the singularities of G(k, lo) in 
the half-plane cj < 0. 

For large k, we see from Eq. (8) that the bounds u>± = (K±k) 2 — 1 are large, and the contribution to negligible, 
except for the lower bound if K is very near k. This makes K also large and, from Eq. (7), one can use the limiting 
form f ~ 1 {K 1 y + uj) ~ A -1 + (ir/2)y/ (K/2) 2 — ui/2, since y is bounded due to the Fermi distribution. The remaining 
integral is easily performed at T = to give: 

T, L (k,w) = -l 1 - (16) 

3 \-i + W2)y/(k/2)*-Q 
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For Lb < this term in the self-energy has no singularity and does not generate any singularity in G(k, u>). Hence we 
would merely have = without the contribution of Sr. 

When T — > the Bosc distribution in Eq. (9) restricts y to negative value, which implies also from the lower bound 
of the integral that K < 2. Hence for large k the logarithm in this equation can be replaced by — 4kK/(k 2 + iu) n ). 
This gives for this case the following expression: 

£ r (/c,u)) = ^L (17) 
u) + k z 

where the constant C is given by: 

r-2 rO 



C = --f dK K 2 [ dy JmT(K,y) (18) 

K JO ip/2-2 



This constant is easily evaluated since the imaginary part of T 1 (K,y) becomes simple from Eq. (7) at T = and, 
if we restrict ourselves to the perturbative situation of small A investigated by Belyakov, its real part is merely A -1 . 
This leads to C = 4A 2 /9. 

Since £r has a pole for u> — — k 2 , it is easy to see that G(k, ui) itself has a pole in the vicinity. In this vicinity Sl 
is small and can be neglected. This leaves us with: 

(T'<k,«) = I* - V - ^ * P + f + (19) 

The pole of G(k, w) at u> — -[k 2 + C'/2k 2 ] has a residue C/(4|^|fc 4 ) which leads to: 

in agreement with Belyakov [22]. One would also obtain this result by performing straight away an expansion in the 
small parameter C, but this is not so obvious to justify since the denominator k 2 + Hi in Eq. (19) can vanish. 



V. DILUTE LIMIT FOR NEGATIVE SCATTERING LENGTH 

The result of Belyakov has been obtained by second order perturbation theory, and the agreement we find with his 
result would let us believe that our approximation is completely valid up to second order in perturbation. This is 
actually not completely correct, as it can be seen by comparing our results with Galitskii's dilute limit theory [23,24]. 
We just sketch here the calculation, at T = 0, which is basically an expansion in powers of the scattering length a of 
the expression of the particle density n in terms of the chemical potential /U. 

If we insert in Eq. (1) the first order expression X^ 1 ) = gn where g = Aira/m, we have for the particle density: 

J (2tt)3 J 27r(iw-&- fl no) 1 ' 

As is well known, this is just the zeroth order result no(n'), with a shifted chemical potential ^ = \i—g n - To proceed, 
we consider the difference n — and we evaluate the self-energy up to second order £(fc, iuj) — gn ~ S^ 2 '(/c, iuj). 
Similarly to Eq. (3) this gives: 



(1) = J _u fL _ i dw £ (2) (M^) 

' ' ' 2tt (iw - - ^ - S( 2 )(fc, iw)) 



I— i 

J (2tt)3 J 



where we have set = — fi' . If one sets directly S^ 2 ' = in the denominator, one finds n — = 0. In order to 
get the proper expansion, we substract from Eq. (22) the corresponding expression without T,^ in the denominator, 
which leads to: 



(1) _ f dk f dcu ^ 2 \k,ito) 
~ J (2tt)3 J 2tt iu- C k 



1 



(23) 



We see that the integrand is important only when iuo — ^ is small, i.e. when the variables are close to uj = and 
k = k' , where = 0. Therefore, one may replace in the integral Y,( 2 \k,ioj) by ^ 2 ^(k' ,0) = C (k' a) 2 fi, where C is 
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a positive constant to be evaluated, and consider a small domain of integration around u> = and k = k' Q . Performing 

the integrations gives n — = — ^(fc a) 2 C, where we have made in this second order term k' w k - We have 

fc 3 

therefore, for n = -^i, the second order expansion we were looking for: 

§2 = ^ ( 2 ™(M - ^o)) 3/2 - (*b «) 2 + • • • (24) 

where k = ^J2m\i and n = k\j(%-K 2 \ 

The constant C" is determined by considering the second order values of the self-energies contributions and 

Ep 2 \ Performing analytically the integrals for Sp \ one finds Sp 2 ^(fc ,0) = (2 + log 2) j^(ko a ) 2 - The integrals 

entering are more easily calculated numerically and we get (fco, 0) = 0.2818 ^z^(koa) 2 . From this we obtain 
the numerical value C" = 0.2597. Inverting the expansion in Eq. (24) in order to express the chemical potential in 
powers of /cjr«, we have: 

k 2 

u = ^{l + ak F a + f3{k F a) 2 + ■■■) (25) 

We get a = ^ and (3 — gfr+C- This is in apparent disagreement with Galitskii's result where (3 = y^^ll— 2 log2) ~ 
0.2597. 

This is due to the fact that our self-consistent calculation docs not include all the second order contributions. Indeed 
when we consider the Hartree term 1771, we see that we have taken it into account only to zeroth order, in Y,^ = griQ. 
However if we want to have a result which is valid up to second order, we want to have the relation n(u) in gn correct 
up to first order while we have it only to zeroth order. Technically this is because we do not have any self-energy 
contribution in the propagator G which comes in the familiar diagram which corresponds to the Hartree term, that is 
we take G = Go in this diagram, instead of taking into account a first order correction to G. We note that this defect 
of our scheme is easily corrected " by hand" : since we know that the exact result for the Hartree term (to any order) 
is anyway gn, proportional to the actual density, we can make directly this correction in our formalism. Indeed if we 
replace in Eq. (24) gno by gn, we obtain (3 = C , in full agreement with Galitskii. 



VI. MOLECULAR BOUND STATE 



As mentionned in section III, we have not yet considered the possibility that T(K, lo) has a pole on the real frequency 
axis, corresponding physically to the existence of a bound state formed by two fermions belonging to two different 
hyperfine states, in other words a molecular state. Naturally this will occur only for a > 0. This is what we will study 
in the present section. 

This pole will appear below the branch point for the cut, which is at K 2 /2 — 2s. Let us call yo{K) the corresponding 
zero occuring in Eq. (7). The equation for Y (K) = yo{K) — K 2 /2 + 2s, which is negative, is given explicitely by: 



r c 

Jo 



dx 



xt/K ln Cosh[((a; + X/2) 2 - S )/2t] 



x 2 -Y /2 cosh[({x- K/2) 2 - s)/2t\ 



(26) 



For increasing coupling strength this zero will first appear when it is right at the branch point, that is lo = 0. This 
gives for the corresponding coupling strength threshold \ t h{K): 



dx 



t l cosh[((x + K/2) 2 - s)/2t] 



xK cosh[({x-K/2) 2 -s)/2t] 



(27) 



As expected this coupling strength is always negative, since one can check that the quantity to be integrated in Eq. 
(27) is positive. This result gives the threshold for formation of a molecule with total momentum K = k$K. In 
particular for zero momentum K = 0, one obtains: 



-Wo)- 



dx 



1 — tanh 



x 2 + l 
2t 



(28) 



which has already been studied in Ref. [15]. In this case we have taken s = — 1 since otherwise one is already in the 
phase diagram domain where the BCS instability occurs [15]. On the other hand for K ^= both s = ±1 are possible. 
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Finally one checks that, for large K Eq. (26) gives merely yo = K 2 /2 — 2s — 8/(7rA) 2 . In particular the threshold 
is given by X^ h (K) = 0. This is expected since in this case atoms become insensitive to the presence of the other 
fermions, and one should recover the result for atoms in vacuum. 

The existence of a bound state implies that we have to consider in Eq. (4) the additional contribution E m coming 
from the corresponding pole occuring at y . It is given by: 



l {k,iu n ) = -i 



dK K fe(^) In 



K 2 



yo 



du 



(29) 



yo 



Introducing: 



J(K) 



f 

Jo 



dx 



cosh[((x + K/2) 2 -s)/2t] 
[x 2 + |F /2|] 2 cosh [((x - K/2) 2 - s)/2t] 



ct/K 



In 



(30) 



we can rewrite the above contribution as: 



kJo 



dK Kb{ 



1 



t ' J(K) 



In 



K 2 



yo 



lU) r , 



(31) 



VII. DILUTE LIMIT FOR POSITIVE SCATTERING LENGTH 



We explore now the behaviour of our approximation on the BEC side of the phase diagram, going to the dilute 
limit. This corresponds to take a > small enough. In this case we will have the presence of molecular bound states, 
corresponding to the pole of the vertex considered in the preceding section. The binding energy becomes large when 
a is small. Since physically the chemical potential will be at most half the bound state energy, it will be large and 
negative in this range. This implies that t = T/\fi\ <C 1 and s — —1. More specifically in this case, the arguments of 
the cosh's in Eq. (26) are always large and we find |Y"o| = 2(2/ nX) 2 . For a molecule at rest (i.e. K — 0) this means 
that the molecular binding energy is H-I^ol = £& = l/(ma 2 ) as expected. Similarly, from Eq. (30), the expression of 
J(K) simplifies into J(K) = n / {4^/\Yo/2) = n 2 \\\/8. 

Turning now to the self-energy we see from Eq. (8) and Eq. (9) that, because of the Fermi and Bose distributions, 
these contributions El and Sr to the self-energy will contain exponentially small factors, respectively exp(— 1/t) and 
exp(— 2/t), since the lower bound for the y integral is larger than (or equal to) 1 or 2, respectively. Hence we can 
neglect these two contributions in this small t limit, and retain only the contribution E m considered in the preceding 
section, which can be checked at the end of our calculation to be indeed exponentially larger than the two other ones. 

The expression Eq. (31) for S m simplifies since the Bose distribution factor forces K to be very small, and the 
argument of the logarithm is very near unity. This leads to: 

W- ^^ (32) 

where the constant: 

C 



is small in the dilute regime we are interested in. 

The expression for the self-energy is made clearer by going back to physical variables and making use of Eq.(36) 
below, which gives: 

v (u ■ n 167ra|^|n/TO(fc o) 2 

Y, m {k,iuj n ) = — ; — — ■ (34) 

«w„ + 4fc + Mb 

where we have set [Ib — 2[i + e^, which we will justify physically just below. Note that T, m (k,iu> n ) has a "hole"-like 
dispersion in the denominator iuo n + ^fe + Ms m contrast to the "particle" -like dispersion iu> n — in the bare Green 
function Go- 

The situation is now similar to the one we had in section IV for the large momentum behaviour. Since we have t <C 1 
(which does not imply that T itself goes to zero), once again only the singularities of G(k,u>) in the half-plane ui < 
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will contribute in the calculation of rik given by Eq. (20). The self-energy Eq. (32) has a pole for u — 1 — |Yo| — fc 2 
and we find that G(k,w) itself has a pole in its vicinity at u) = 1 — |Yo| — fc 2 — C/(|Yo| + 2fc 2 ). The corresponding 
residue is C/(\Y \ + 2k 2 ) 2 . This leads to: 

nk = (\Y \ + 2fc 2 ) 2 (35) 
and after integration over k we have finally for the total density for a single hyperfine state: 



n _ C r 



k 2 C / 2 C 

dl m+ 2k^ = wJwr^ koa (36) 



These results Eq. (35) and Eq. (36) are just what one expects in this regime. Indeed taking Eq. (33) into account 
and going back from reduced to physical variables we can rewrite Eq. (36) as: 

47r f°° k 2 

n = j- ,„ / dk t> (37) 

(27t)3 7 e xp[(£ _ eb 2m )/T] - 1 

which is just what is expected for the density of non interacting bosonic molecules with binding energy ej,, mass 2m 
and chemical potential \ib = 2[i + ej. Similarly we know that the wavefunction for the relative motion of atoms in 
these large nearly resonating molecules is ip(r) = (A/r)e~ r l a , where A is a normalization constant. The momentum 
distribution is proportional to the square of the Fourier transform of this wavefunction, that is (fc 2 -I- l/a 2 )~ 2 , which 
is just what we find in Eq. (35) in reduced variables, since \Y$\ = e&/|/x| = 2/(fcoa) 2 . It is worthwhile to note that we 
are in the normal state, in contrast with the similar result obtained from the approximate BCS wave function in the 
strong coupling limit, which holds in the superfluid state. 

We may wonder if our calculation could also produce interactions between our molecular states. Although this is 
not so obvious, this does not seem to be the case. Indeed in our expression Eq. (31) for the self-energy, we have only 
a single Bosc distribution which is appearing, which merely leads to the expression Eq. (36) for the density where 
this Bosc factor appears again. Hence if we had interactions we would expect the appearance of products of Bosc 
distributions, which we do not have. In a related way we notice that this Bose factor is linked to the existence of 
ladders in the diagrammatic writing of our formulation. This corresponds physically to the propagation of a single 
molecule. However in our formulation there are no diagrams where two ladders interact. Physically this means that 
we have no interactions between molecules. In order to find them one has to consider more complicated diagrams, 
describing such an interaction, as it has been done for example by Pieri and Strinati [16] and very recently by Brodsky 
et al [27] , giving an exact agreement with the four fermions calculation of Petrov et al [28] . 



VIII. TEMPERATURE EVOLUTION OF THE SYSTEM IN THE MOLECULAR LIMIT 

In this dilute molecular limit £b = l/(ma 2 ) 3> E F it is interesting to consider in more details the temperature 
evolution of the system. At high temperatures T 3> E F the situation is governed here by the dynamical equilibrium 
between molecules and unbound fermions. This situation is described by the well-known Saha (or law of mass action) 
formula (see [20]). In 3D it reads: 

li= * ( m T) 3 / 2 exp(-e b /T), (38) 
ub ^tt 

where the total particle density n = 2n = n F + 2n B , with n F the sum of the free fermion density and ns the bosonic 
molecular density. The crossover temperature T* for which n F = 2n b = n/2 is given by: 

E (39) 



2 111 E F 



where the logarithm in denominator of Eq.(39) has an entropic character. For T » T, one has n F S> 2ub, hence 
nw n F and the fermionic chemical potential \i m —^T\n{T /T^ EC ) has the standard Boltzman form (T C BBC ~ E F is 
a typical temperature for Bose-Einstein condensation). 

At much lower temperatures E F <C T <C T* the situation drastically changes. The number of unpaired fermions 
n F ~ exp(— eb/2T) becomes exponentially small and hence h ~ 2ns- The fermionic chemical potential acquires a 
kink: 
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^- e i-\ Tln T^c (4°) 

In this low temperature regime we have essentially 2[i = — eb, that is koa = 1. This simplifies the expression Eq.(34) 
into: 

„ . I6ira\u\n/m , . 

S ro (k, iu n ) = -Vp-L 41 

«W„ + i k + (J, B 

It is important to emphasize that for T<T, our two-particle vertex has the simple pole structure: 

T(K, iw v ) = , 4| ^r" , (42) 

where, as we have seen in the preceding section, /xg = 2/j + is the molecular chemical potential. Correspondingly 
we have \ib ~ — §Tln ■ 

The dressed one-particle Green function G _1 (k, = G^ 1 (k, iuj n )~ S(k, iw n ) has a two-pole structure for T <C £?f: 

G = \ fi | | , , (43) 

where, as before, = k 2 /2m — /j,. 

The spectral function A(k, uS) = — — ImG(k, u> + i0+) reads: 

</, \ ,„ 47ra|u|n/m. , . A-Ka\u\n/m . , , 

4(k, w) « (1 !g_i_)5( w - a) + 75-^^ + & + Ms) (44) 

It reflects for r«T, the existence of two bands: the filled bosonic band and, separated by the correlation gap ef,, 
the almost empty band of unbound fermions. Integrating the spectral weight it is easy to check, that in this regime: 



4Tra\fi\n 1 f°° k 2 dk 



J 

Jo 



m 2n 2 J Q e k 

The specific heat of the system: 

dE d , f k 2 k 2 dk , k 



— n w 2n B (45) 



dE d . r k 2 k 2 dk , k 2 N /i Bl 

= t^; = t^; / -; -^-cxp -- — — I exp — - ~ n = const. (46) 

dT dT l J Am 2tt 2 i y AmT 1 T l v ' 

is temperature independent in agreement with general thermodynamic requirements. 

IX. GENERAL MOMENTUM DISTRIBUTION 

An intermediate step in our calculation is naturally the particle momentum distribution n k = n$k + Snk, where 
nok = (cxp(£fc/T) + 1) _1 is the Fermi distribution and the correction due to interactions is, with our reduced units, 
given by: 

k ^ S [i^n - k 2 + s - Z(k,iw„)].[iu n - k 2 + s] ^ ^ 

with k = kok and we have used — iu) n ) = S*(fc,iw n ), and we have explicitcly in terms of the contributions 
considered above Y,(k,iu) n ) — T,r(k,iu n ) + S_L(fc,i(D„) + Y, m (k,iui n ). Then the particle number is obtained by n = 
(1/27T 2 ) J^°dk k 2 Ufc. The calculations of Srik and of n are both handled numerically. 

The above Matsubara summation is numerically quite convenient since it converges fairly rapidly (for large ui n the 
terms in the series Eq. (47) behave typically as 1/w^). Nevertheless it is of interest to consider another possible 
calculation where this summation would be transformed in a contour integration over frequency as we have done for 
example in section IV and VII. The contour would be deformed to enclose the singularities of the Green's function 
occuring on the real frequency axis, and the result could be expressed in terms of the spectral density. This would 
allow to ascribe a physical meaning to the various contributions, as we have done in section IV where we had to deal 
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with a fermion pole, whereas in section VII we had a molecular pole. This is of interest since the total particle number 
n is very often split in the literature into a free fermion term plus a molecular term asn = nf erm + n mo i . 

It is easy to see that this split is not possible in general because the Green's function has a cut on the real axis 
extending from — oo to +00. There is no way to split the Green's function itself in a sum of a free fermion term and a 
molecular term since this is rather at the level of the self-energy that such a separation occurs. We have indeed found 
for the self-energy a molecular contribution S TO , linked to the molecular pole of r(k, u>). The contribution Sr can also 
be understood as linked to molecules (as it is clear from the Bose factor it contains in Eq. (9)), but this is rather the 
continuous spectrum of molecules broken into two fermions, rather than the discrete spectrum linked to S TO . Finally 
the contribution is linked to single fermions since it arises in Eq.(4) from the pole of the Green's function Go 
which is clearly signaled by the Fermi factor in Eq.(8). We could still think of separating a molecular contribution 
from a fermionic contribution if we had in the spectral density of the Green's function an energetically well separated 
part, similar to what we had in section VII. However this is likely to occur only when the molecular binding energy 
is large compared to temperature, in which case we will have dominantly molecules, and the separation of the density 
into a molecular part and a fermionic part is rather uninteresting. Anyhow it is clear from the above discussion that 
such a separation is in general unwarranted and, if it is used, it must taken cautiously. 



X. COMPARISON WITH OTHER WORKS 



As we have already indicated our approach is in complete agreement with the equations written by Perali, Pieri, 
Pisani and Strinati [17] when we restrict them to the normal state. It is now also of interest to compare our framework 
with the one used by Nozieres and Schmitt-Rink [7] (NSR). 

Quite generally we can write [25] the thermodynamic potential (function of \i and T) of the two hypcrfine states to 
as: 

/■A j\i 

n-n Q = — T y~] Gq 1 (k, uj n )[G(k, uj n ) — G (k, u n )\ (48) 

where the integration is over a varying coupling constant A'. Here CIq is the thermodynamic potential for the non 
interacting gas, G(k, w n ) and Go(k, uj n ) arc the Green's functions for the interacting and the non interacting gas. 
Introducing the expansion of the full Green's function in terms of the self-energy G = GoE^ot^^)''' ^ s can ^ c 
rewritten as: 

n-n = I ^rT V V[s(k, w „)G (k,c„)] p (49) 



2 °=/ ^r^]T[s(k,^)Go(k,c„)r 
J ° k, WnP =i 



In our approximate scheme the self-energy is given by Eq.(4), with T = — 2ir 2 T /(mko) and f _1 = A -1 + / from Eq. 
(7), where / denotes the integral in the right-hand side of this equation. If we consider the first order term SGo in 
the above expansion we can rewrite its contribution as: 

T^S( kiUn ) Go (k, Wll ) = r^ iTW^M E(- 1 ) P ~ 1 [A/(k>;)F (50) 

where w' n is a bosonic Matsubara frequency. To obtain this equation we have used T w ^ Go(k, w„)Go(k' — k, u)' n — 
uj n ) = —mk n I(k',u>' n )/(2Tr 2 ), as it results from comparing Eq.(6) and Eq.(7), and we have omitted the regularization 
term of the high momentum divergence, since NSR use rather a separable potential to achieve this regularization. 
Inserting this contribution in Eq.(49) we find for the corresponding contribution to the thermodynamic potential: 

00 (— Fl p_1 

n - Oo = t E — — [ A/ ( k < = T E M 1 + A/ ( k > ( 51 ) 

k,oi n p=l k,w n 

which is just NSR expression, provided that we identify their x(k, u> n ) with our — A/(k, u) n ), as it can be confirmed 
by considering their detailed expressions. Accordingly NSR calculations take into acount only the first term in the 
expansion Eq.(49), and uses the same approximate self-energy as we do. 

If one calculates directly the particle number by Eq.(l), instead of using the thermodynamic potential f2 and 
n = —dfl/d/j,, the NSR scheme corresponds to take G — Go = GqSGo, that is again the first order term in the above 
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expansion of G, as it can be seen by performing the differentiation with respect to fi on Eq.(51). By contrast our 
expression for G in Eq.(l) contains all orders in this expansion of G. Hence our scheme goes one step further by 
retaining all the terms in the expansion, while naturally keeping the same approximate expression for the self-energy. 
In other words we allow any number of self-energy inclusions in our propagator. We note that, in the dilute limit, it 
is perfectly valid to retain only the first order term, as NSR do. This is coherent with the fact that they find properly 
the zeroth order expression for the Bose condensation temperature. Diagrammatically they have for the propagator 
the familiar single "wheel" contribution (i.e. a ladder diagram with one side closed on itself). In our case we have 
any number of wheels interconnected as can be seen in Fig. 1. Physically this means that, in our case, the molecular 
bosons corresponding to the ladder diagrams can be broken and the elementary fermions present in the molecule can 
reappear individually. 





XI. CRITICAL TEMPERATURES AND MOLECULAR INSTABILITY 

In the case of the BCS transition temperature, for which yU > 0, the relation between the critical temperature T c 
and the chemical potential is given by the standard equation: 



1 

A 



dx 



o 



x 2 

1 — — ^ tanh ■ 



2t c 

where t c = T c / fj,. With the standard physical variables this equation reads: 



(52) 



O poo 

l = f£/ dk [l-^^tanh^^] (53) 
7T J 1 e k - n 2T C J V 1 

where the chemical potential [i is positive. This equation is obtained by writing that T(K, uj) has a pole at the chemical 
potential for zero total momentum K = 0, that is by setting Yq = 2, s = 1 and K — ► in Eq. (26). This BCS transition 
line terminates when fi = (this implies t c — > oo) which gives the relation between the limiting critical temperature T 
and the corresponding scatt ering leng th a . This is more directly obtained from the equation with physical variables 
Eq.(53) and is given by tt / y/^ma^To — J °° dz[l — tanh(2 2 /2)], i.e. explicitely ma^To = 1.07 = T /ebo in terms of the 
molecular binding energy e^o = I/mag at this point. Obtaining T /Ep and 1/kpao requires the numerical calculation 
of the gas density at that point (knowing that it satisfies [i = 0), along the procedure described in the preceding 
section IX. 

In the weak coupling limit \a\kp -c 1 this gives for the dominant order T c <~ Ep exp(— 1/A). When the calculation 
is carried out to the next order term, one finds the numerical coefficient in front of the exponential as: 

T-^0.61^ex P (-^) (54) 

However, already at the level of this coefficient, the simple weak coupling limit is not correct and one needs to take 
into account polarization contributions to the effective interaction, arising from second order diagrams in the gas 
parameter A, as it has been shown by Gor'kov and Mclik-Barkhudarov [30]. This leads to their well-known formula: 
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T-^0.28^exp(-^) (55) 

Hence, at this level already, our approximation does not yield the exact result. In order to improve it on this side 
we should take into account polarization diagrams at the level considered by Gor'kov and Melik-Barkhudarov, and 
beyond as it has been done for example in Ref. [31], [32] and [33]. 

On the other hand we meet also naturally in our calculations the Bose-Einstein transition for the molecular states 
we have discussed in sections VI and VII. Indeed it is reached when the Bose factor b(yo/t) entering the self-energy in 
Eq. (31) diverges, since this factor corresponds physically to the statistical occupation of the molecular states. This 
divergence occurs for yo{K) — 0. Naturally this is for a zero total molecular momentum K — that it occurs first. 
This implies that Yq(K) = yo(K) — K 2 /2 + 2s = —2, since we need s = — 1 in order to have K = molecules, as 
noted above. From Eq. (26) we find for the critical temperature T c = \fj,\t c : 



I />oo 

dx 



|A| 



o 



x 2 x 2 
1 — — tanh ■ 



1 2t r _ 



(56) 



where we use the fact that this transition occurs only for A < to write —A = |A|. We see from this equation that, 
by letting n — ► (which implies again t c — > oo), the Bose-Einstein transition line terminates at the same point [i = 
and T — To as the BCS transition line. This Eq.(56) is actually in agreement with the finding of Ref. [8] (as more 
easily seen when we go back to unreduced units). However, in contrast to the approach of these authors which treat 
the superfluid state, our reasoning involves only the normal state and its instabilities. 
Going back to physical variables the above equation reads: 

1 = * f dk {1 _ tanh f* + N ] (57) 

Wo 1 e k + 2T C J V 1 

In the dilute molecular limit kpa <C 1 (or equivalently na 3 <C 1), we have 3> Ep and we have seen in section VII 
that fiB = 2[i + €b, implying \fj,\ > e(,/2. This gives in Eq.(57) \fi\ 3> T c ~ Ep and hence tanh[(efe + |/i|)/2T c ] w 1. 
Accordingly: 



TT JO e k + |Ml 



(58) 



Hence in this limit our equation reduces to Hb(T c ) = as expected for the Bose-Einstein condensation of dilute gas. 
Using the standard formula for bosons [20] with density ns — n and mass 2m, one finds [29]: 

Tf EC = 3.31 ^ 0.218E F , (59) 

Actually in this limit, our approximation is somewhat deficient as we have noted above in section VII. In fact 
in this regime T < T, we are dealing with a low-density bose-gas of molecules with effective hard-core repulsion. 
The molecule-molecule scattering length clm — 0.6 a has been first calculated recently by Petrov et al [28] from the 
4-particle Schrodinger equation and recently rederived diagrammatically by Brodsky et al [27]. The inclusion of 
molecule-molecule interaction leads in this case to a correction to critical temperature: 

T C BEC = 0.218 E F [1 + a(a M n 1/3 )] (60) 

where a ~ 1.3 is a numerical coefficient obtained by Kashurnikov, Prokof'ev, and Svistunov [34] in the framework of 
rather complex Monte-Carlo calculations. 

It is also worthwhile to note that actually, in the general case, we do not need the above physical interpretation 
of the Bose factor to find the Bose-Einstein transition. Indeed if we would calculate as discussed above the particle 
number as a function of temperature and chemical potential, and if we would enter the condensate domain of the 
phase diagram without realizing it, we would find that we could not accomodate all the particles we have in our 
system, even by further lowering the temperature, because the particle number would decrease with temperature. 
This would signal the appearance of a condensate in order to accomodate the remaining particles. In other words we 
can not miss the transition, which appears automatically in our approach. 

Finally our major interest is the line in the phase diagram where molecules begin to form. As it has been stressed 
in Ref. [15] this threshold line is shifted from its standard location for molecules in vacuum, which is at unitarity. This 
is due to the presence of all the other fermions, which influence the formation of a molecule in much the same way as 
they do for the formation of Cooper pairs. Moreover, in contrast with the situation in vacuum and in similarity with 



14 



Cooper pairs, the location of the threshold line where molecules begin to form depends on the total momentum of 
the molecule under consideration. For a very large momentum the molecule is insensitive to the presence of the other 
fermions, and the threshold line is at unitarity a" 1 = 0, just as in vacuum. On the other hand the molecules with 
zero momentum are the most sensitive to the presence of the other fermions. The question of molecular formation 
has already been considered in section VI and the location of this threshold line is given by Eq.(27). Again this line 
terminates at /i = 0, at the same point as the BCS line and the BEC line. Hence the three physical lines of interest 
in this section terminates at the same point. 




-1 1 

1 /k p a 

FIG. 2. Threshold line T mo i for the formation of molecules with zero total momentum K = 0, critical temperature for the 
BCS instability T(, cs and critical temperature for the Bose-Einstein condensation of molecules Tf, ec , as a function of the scattering 
length a. Dots are the actual numerical results. Curves are smooth interpolations through these results. The wavevector Uf is 
related to the gas density n by n = k F /(6iT 2 ), and the Fermi energy is defined by Ef = h 2 k F /(2m). The line ft = is drawn 
to show roughly the region (below this line) where the physics is qualitatively the one of a degenerate Fermi gas with pi > 0. 
The terminal point ao,To common to the molecular line, the BCS line and the BEC line is located at l/(kFdo) = 0.37 and 
T q /Ef = 0.29. The BEC limit for the dilute molecular gas is at T c = 0.218£ F 



The last step in order to draw a physical phase diagram is to express the critical temperature in terms of the 
particle density n, rather than in terms of the chemical potential \i. This is done by calculating numerically n as a 
function of fi, as indicated in section IX. As indicated in the introduction, we use as usual wavevector and energy 
scales directly related to the atomic density n for a single hyperfine state, namely the Fermi wavevector kp defined 
by n = fc^/(67r 2 ), and the Fermi energy Ep = k F /(2m). The results are displayed in Fig. 2. 

Let us first discuss the BCS and the BEC lines. Qualitatively the general shape is similar to the NSR result. 
The maximum of T c /Ep is very close to the terminal point ao,T and (at the scale of our figure) it is very smooth. 
This shape is also found in later work [8,17], including phcnomenological approaches [29]. Quantitatively our results 
are slightly higher than those of Perali et al [17] for the homogeneous case, as it can be seen from the in their 
Fig.l. Although this discrepancy is rather small, it is somewhat puzzling since both numerical calculations have been 
performed in a careful way, although the specific starting equations for the numerics are different as we have seen in 
section II. Anyway the difference can be taken as a typical uncertainty arising in this kind of numerical calculations. 

We consider finally the molecular threshold line. The quite interesting feature is that, for all the temperatures 
satisfying T < Ep (and clearly even somewhat above), 1/kpa stays typically around 0.4, that is roughly its value at 
the terminal point. This is in contrast with what one would expect from the fact that, at high temperature, this line 
goes toward the unitarity line 1 /kpa = since this quantum shift of the molecular threshold clearly disappears in the 
classical regime. Hence the region between the molecular line and unitarity could have been pretty small, in contrast 
with what we find in Fig. 2. This makes naturally easier the experimental observation of the shift. It may even be that 
it has been already seen in the very recent vortex experiment [5]. Indeed in order to observe vortices, experimentalists 
have to let the Fermi gas expand while at the same time forming molecules, in order to see the depression in molecular 
density associated with the vortex in the same way as it is done for standard Bose condensates. In the experimental 
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process, it has been observed that, when forming vortices on the BCS side a < 0, it is necessary to switch the magnetic 
field to somewhere below the Feshbach resonance in order to be able to see the vortices (the corresponding l/kpi 
is of order 0.35, in nice agreement with the location of the molecular line). A possible explanation might be that, 
because (low momentum) molecules can not form beyond the molecular line, the formation of molecules necessary 
to see the vortices can not occur when the magnetic field is not brought to low enough values. However the answer 
lies in the trajectory of the Fermi gas in Fig. 2 during expansion, which is not at all obvious to determine since the 
expansion is clearly a complicated dynamical process. Hence we can certainly not exclude that the explanation for 
this experimental observation lies somewhere else. 



XII. CONCLUSION 



In this paper we have treated in a self-consistent way the various molecular instabilities arising in the normal state 
of a degenerate Fermi gas. We have stressed that these instabilities have to be contrasted with the smooth crossover 
in the nondegenerate regime, corresponding to the temperature where the formation of the predominant molecules 
occurs as the temperature is lowered. We have covered all the range of scattering lengths, so as to cover the whole 
BEC-BCS crossover. The molecular instabilities manifest themselves mostly as poles of the vertex corresponding 
to particle-particle scattering. This vertex has been calculated in a ladder approximation. The threshold for the 
appearance of standard molecules corresponds to the existence of a pole at zero energy in the vertex. The BCS 
instability corresponds also to the appearance of molecular-like objects, namely Cooper pairs, but the pole in the 
vertex appears at the chemical potential. Finally the equation giving the Bose-Einstein condensation instability is a 
simple, but elegant, continuation of the equation giving the BCS instability. In order to find the critical temperatures 
corresponding to these various instabilities, we have taken into account the interactions between fermions in this 
normal gas, by making use of the same vertex as the one used to obtain the instabilities themselves. This leads to the 
self-consistent T-matrix approximation. We have shown that this approximation is quite satisfactory since it leads 
in a number of limiting cases to results which are in agreement with known exact results. These are specifically the 
high temperature regime, the large momentum limit at T = 0, the dilute limit both on the BCS and the BEC side of 
the crossover, which we have investigated successively. The calculated phase diagram shows that, in particular, the 
threshold for formation of molecules at rest undergoes a sizeable shift toward the BEC side, due to the hindering effect 
of the quantum gas on the molecular wave function. This shift remains important up to temperatures comparable to 
the Fermi energy of the gas. Finally our approximation is somewhat defective on the BEC side since it does not allow 
to describe interaction between molecules. Improvement is clearly needed in this direction and will be the subject of 
separate publication. 
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